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Abstract. Wc consider homogenization for weakly coupled systems of Hamilton-Jacobi 
equations with fast switching rates. The fast switching rate terms force the solutions 
converge to the same limit, which is a solution of the effective equation. We discover the 
appearance of the initial layers, which appear naturally when wc consider the systems 
with different initial data and analyze them rigorously. In particular, we obtain matched 
asymptotic solutions of the systems and rate of convergence. We also investigate properties 
of the effective Hamiltonian of weakly coupled systems and show some examples which do 
not appear in the context of single equations. 
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1. Introduction 

In this paper we study the behavior, as e(> 0) tends to 0, of the viscosity solutions 
[u\ , uf) of the following weakly coupled systems of Hamilton- Jacobi equations 



( (u\) t + Hi(~, Du\) + j(ul -u%) = in R n x (0,T), 
(u|) t + tf 2 (~ Dul) + j(u e 2 - u\) = in W 1 x (0, T), 



k uf(z, 0) = fi(x) on M n for i = 1,2, 

where T > 0, ci, c 2 are given positive constants and the Hamiltonians Hi(£,p) : lR n x M" — > 
R are given continuous functions for % = 1,2, which are assumed throughout the paper to 
satisfy the followings. 

(Al) The functions Hi are uniformly coercive in the ^-variable, i.e., 

lim inf{#i(£,p) ^ G R", \p\ >r} = oo. 

r—toc 

(A2) The functions £ H- Hi(£,p) are T n -periodic, i.e., Hi(£ + z,p) = Hi(£,p) for any 
£,p el",26 Z n and i = 1,2. 

The functions j] are given continuously differentiable functions on MJ 1 with ||-D/j||L°°(Rn) 
are bounded for i = 1,2, respectively. Here u\ are the real- valued unknown functions on 
]R n x [0,T] and {uf) t := duf/dt, Duf := (duf/dxi,...,duf/dx n ) for z = 1,2, respectively. 
We are dealing only with viscosity solutions of Hamilton-Jacobi equations in this paper 
and thus the term "viscosity" may be omitted henceforth. 

1.1. Background: Randomly Switching Cost Problems. System (C e ) arises as the 
dynamic programming for the optimal control of the system whose states are governed by 
certain ODEs, subject to random changes in the dynamics: the system randomly switches 
at a fast rate 1/e among the two states. See [13, 14, 16] for instance. Also see [26, 7] for 
another switching cost problems. In order to explain the background more precisely, we 
assume in addition that the Hamiltonians Hi are convex in p here. We define the functions 
u\ : W 1 x [0,T] ->■ R by 

ut(x,t) :=inf{E,(^ ^(.)(^, -*)(*)) da + U(t)(v(t))) }, (1.1) 

where L, : M. 2n — > M U {+00} are the Fenchel-Legendre transform of ifj, i.e., Lj(£,g) := 
su PpeK™ G° ' <? — Hi(£,p)) for all (£,5) G M 2n and the infimum is taken over 77 G AC ([0,t],R n ) 
such that 77(0) = x. Here AC ([0,t],M n ) denotes the set of absolutely continuous functions 
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with value in M. n and Ej denotes the expectation of a process with z/ £ (0) = % where v is a 
{l,2}-valued process which is a continuous-time Markov chain such that 

P(z/(s + As) = j | z/(s) = i) = jAs + o(As) as As -> for i ^ j, (1.2) 

where o : [0, oo) — > [0, oo) is a function which satisfies o(r)/r -)• as r -)• 0. Formula 
(1.1) is basically the optimal control formula for the solution of (C e ), where the random 
switchings among the two states are governed by (1.2). 

We first give a formal proof that (u\ , w|) given by (1.1) is a solution of (C e ). The rigorous 
derivation will be proved in Appendix by using the dynamic programming principle. We 
suppose that Ui G C 1 (lR n x [0, T]) here. Set u(x,i,t) := Ui(x,t) and Y(s) := (r/(s), z/ £ (s)) 
for t] G AC (R n ) with r/(0) = x and let v £ be a Markov chain given by (1.2) with i> £ (0) = i. 
By Ito's formula for a jump process we have 

E i (u £ (Y(t),0)-u £ (Y(0),t)) 
= E i(J o - u l( Y ( s )^ ~ s ) + Du £ {Y{s),t - s) • r]{s) ds 

+ f E (« e fo( a U *) - «'fo00, ^00, s)) ■ ^ ds) 

J ° 3=1 £ ' 

>Ei( I -ul(Y(s),t-s)-H u . (s) ( 1 l,Du £ )-L uHs) ( 1 l,-fi)ds 
^ Jo £ £ 

+ f E ("'M'U s ) - uE (v(s), S(s), s)) • ^ ds) 

y ° i=l 5 J 

= -Ei(^J^ L v . {s) (^-,-7])ds^j. 

Thus, 

u £ {x,i,t) < E *(^ L u e {s) {^-r])ds + u £ (Y(t),0)y 

In the above inequality, the equality holds if — r)(s) G D~H v e^{r]{s)/e ) Du(Y(s), s)), 
where D~Hi denotes the subdifferential of Hi with respect to the p-variable. 

1.2. Main Results. There have been extensively many important results on the study 
of homogenization of Hamilton-Jacobi equations. The first general result is due to Lions, 
Papanicolaou, and Varadhan [27] who studied the cell problems together with the effective 
Hamiltonian and established homogenization results under quite general assumptions on 
Hamiltonians in the periodic setting. The next major contributions to the subject are due 
to Evans [16, 17] who introduced the perturbed test function method in the framework 
of viscosity solutions. The method then has been adapted to study so many different 



4 



H. MITAKE AND H. V. TRAN 



homogenization problems that we cannot describe the complete list of references. We here 
only refer to the papers related to our work. Concordel [11, 12] achieved some first general 
results on the properties of the effective Hamiltonian concerning flat parts and non-flat 
parts. Afterwards Capuzzo-Dolceta and Ishii [9] combined the perturbed test functions 
with doubling variables methods to obtain the first results on the rate of convergence of 
u £ to u. We refer to [5, 35] for some recent progress. 

There have been some interesting results [33, 6, 8] on the study of homogenization for 
weakly coupled systems of Hamilton-Jacobi equations in the periodic settings or in the al- 
most periodic settings. We refer the readers to [15, 23] for the complete theory of viscosity 
solutions for weakly coupled systems of Hamilton-Jacobi and Hamilton-Jacobi-Bellman 
equations. Since the maximum principle and comparison principle still hold, homogeniza- 
tion results can be obtained by using the perturbed test function method quite straightfor- 
wardly with some modifications. Let us call attention also to the new interesting direction 
on the large time behavior of weakly coupled sytems of Hamilton-Jacobi equations, which 
is related to homogenization through the cell problems. The authors [30], and Camilli, 
Ley, Loreti and Nguyen [7] obtained large time behavior results for some special cases but 
general cases still remain open. 

Let us also refer to one of the main research directions in the study of homogenization, 
stochastic homogenization of Hamilton-Jacobi equations, which were first obtained by 
Souganidis [34], and Rezakhanlou and Tarver [31] independently. See [28, 25, 32, 29, 1] for 
more recent progress on the subject. 

First we heuristically derive the behavior of solutions of (C e ) as e tends to 0. For 
simplicity, from now on, we always assume that C\ = c 2 = 1. We consider the formal 
asymptotic expansions of solutions of (C e ) of the form 



Set £ := x/e. Plugging this into (C e ) and performing formal calculations, we achieve 
(ui) t + . . . + Hi(£, D x Ui + H ) + -(ui-Uj) + (vi — Vj) H = 0, 

where we take i,j G {1,2} such that {i,j} = {1,2}. The above expansion implies that 
u\ = u 2 —'■ u. Furthermore, if we let P = Du(x,t) then (vi,v 2 ) is a T n -periodic solution 
of the following cell problem 



where H(P) is a unknown constant. Because of the T n -periodicity of the Hamiltonians 
Hi, we can also consider the above cell problem on the torus T n , which is equivalent to 



u\{x, t) := Ui(x, t) + evi(-) + 0{e 2 ). 



(E P ) 



H x {t, P + D Vl (C, P)) + P) - V2& P) = H(P) 
# 2 (£, P + Dv 2 (£, P)) + v 2 (£, P) - vi(£, P) = H(P) 



in M n , 
in M n , 
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consider it on M. n with T n -periodic solutions. By an argument similar to the classical one 
in [27], we have 

Proposition 1.1 (Cell Problems). For any P e W 1 , there exists a unique constant H(P) 
such that (Ep) admits a T n -periodic solution (i>i(-, P), v 2 (-, P)) E C(lR n ) 2 . We call H the 
effective Hamiltonian assosiated with {Hi^Hz). 

See also [4, 30, 7] for more details about the cell problems for weakly coupled systems. 
Our main goal in this paper is threefold. First of all, we want to demonstrate that u\ 
converge locally uniformly to the same limit u in IR n x (0, T) for i — 1, 2 and u solves 



This part is a rather standard part in the study of homogenization of Hamilton- Jacobi 
equations by using the perturbed test function method introduced by Evans [16] with 
some modifications. The only hard part comes from the fact that we do not have uniform 
bounds on the gradients of u\ here because of the fast switching terms. We overcome this 
difficulty by introducing the barrier functions (see Lemma 2.1) and using the half- relaxed 
limits (see the proof of Theorem 1.2). The barrier functions furthermore give us the correct 
initial data for the limit u. Let (wf,wf) be the solution of (C e ) henceforth. 

Theorem 1.2 (Homogenization Result). The functions u\ converge locally uniformly in 
R n x (0, T) to the same limit u G C 0,1 (K n x [0, T]) as e ->• for i = 1, 2 and u solves 



Formula (1.1) of solutions of (C e ) actually gives us an intuitive explanation about the 
effective initial datum /. As we send e to 0, the switching rate becomes very fast and 
processes have to jump randomly very quickly between the two states with equal probability 
as given by (1.2). (Note that we are assuming c\ — c 2 — 1 now.) Therefore, it is relatively 
clear that / is the average of the given initial data /j for % — 1, 2. 

The second main part of this paper is the study of the initial layers appearing naturally 
in the problem as the initial data of u\ and u are different in general. We first study 
the initial layers in a heuristic mode by finding inner and outer solutions, and using the 
matching asymptotic expansion method to identify matched solutions (see Section 3.1). 
We then combine the techniques of the matching asymptotic expansion method and of 
Capuzzo-Dolceta and Ishii [9] to obtain rigorously the rate of convergence result. 

Theorem 1.3 (Rate of Convergence to Matched Solutions). For each T > 0, there exists 



u t + H(Du) =0 inf" x (0, T). 




(1.3) 



C := C{T) > such that 



ml 



£°°(K n x[0,T]) 



< Ce 1/3 fori = 1,2, 
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where u is the solution of (1.3) and 

ml(x, t) := u(x, t) + (/ *~ 2 /j)(x) e-y (1.4) 
with j G {1,2} such that {i,j} = {1,2}. 

Finally, we study various properties of the effective Hamiltonian H. It is always ex- 
tremely hard to understand properties of the effective Hamiltonians even for single equa- 
tions. Lions, Papanicolaou and Vadrahan [27] studied some preliminary properties of the 
effective Hamiltonians and pointed out a 1-dimensional example that H can be computed 
explicitly. After that, Concordel [11, 12] discovered some very interesting results related 
flat parts and non-flat parts of H for more general cases. Evans and Gomes [18] found 
some further properties on the strict convexity of H by using the weak KAM theory. 

The properties of H for weakly coupled systems of Hamilton-Jacobi equations in this 
paper are even more complicated. In case Hi = if 2, the effective Hamiltonian for the 
weakly coupled systems and the single equations are obviously same. Therefore, we can 
view the cases of single equations as special cases of the weakly coupled systems. However, 
in general, we cannot expect the effective Hamiltonians for weakly coupled systems to have 
similar properties like single equations' cases. 

The first few results on flat parts and non-flat parts of H are generalizations to the ones 
discovered by Concordel [11, 12], and are proved by using different techniques, namely the 
min-max formulas which are derived in Section 4.2 and the constructions of appropriate 
subsolutions. On the other hand, we investigate other cases which show that the properties 
of the effective Hamiltonians for weakly coupled systems are widely different from those 
of the effective Hamiltonians for single equations. Theorems 4.14, 4.17, 4.18 , 4.20, which 
are ones of our main results, describe some rather new results which do not appear in the 
context of single Hamilton-Jacobi equations. Since the theorems are technical, we refer the 
readers to Section 4.3 for details. 

We are grateful to L. C. Evans for his suggestion which leads us to this project. We 
thank G. Barles, D. Gomes, H. Ishii, T. Mikami, and F. Rezakhanlou for their fruitful 
discussions. We also thank S. Armstrong and P. E. Souganidis for letting us know about 
the coming result on stochastic homogenization of weakly coupled systems of Hamilton- 
Jacobi equations of B. Fehrman [19]. Fehrman [19] independently obtained interesting 
homogenization results of monotone systems of viscous Hamilton-Jacobi equations, which 
are similar to ours, with convex Hamiltonians in the stationary, ergodic setting by using 
the ideas of Armstrong and Souganidis [1, 2]. His work includes as well generalizations to 
other related systems. 

The paper is organized as follows. In Section 2 we prove the homogenization result, 
Theorem 1.2. Section 3 devotes to the study of initial layers and rate of convergence. We 
derive inner solutions, outer solutions, and matched asymptotic solutions in a heuristic 



HOMOGENIZATION OF WEAKLY COUPLED SYSTEMS OF HJ EQUATIONS 



7 



mode and then prove Theorem 1.3. The properties of the effective Hamiltonian are stud- 
ied in Section 4. We obtain its elementary properties in Section 4.1, the representation 
formulas in Section 4.2, and flat parts, non-flat parts near the origin in Section 4.3. In 
Section 5 we prove generalization results for systems of m equations for m > 2. We then 
prove also the homogenization result for Dirichlet problems and describe the differences 
of the effective data between Cauchy problems and Dirichlet problems in Section 6. Some 
lemmata concerning verifications of optimal control formulas for the Cauchy and Dirichlet 
problems are recorded in Appendix. 

Notations. For k G N and A C M n , we denote by C(A), C 0,1 (A) and C k (A) the space of 
real- valued continuous, Lipschitz continuous and fc-th continuous differentiable functions on 
A, respectively. We denote L°°(A) by the set of bounded measurable functions and ||-||l<x>m) 
denotes the superemum norm. Let T n denote the n-dimensional torus and we identify T n 
with [0, l] n . Define II : ~R n — > T n as the canonical projection. By abuse of notations, we 
denote the periodic extensions of any set BcT" and any function / G C(T n ) to the whole 
space M. n by B, and / themselves respectively. For a, b G M, we write aAb = min{a, b} and 
a V b = max{a, b}. We call a function m : [0, oo) — > [0, oo) a modulus if it is continuous, 
nondecreasing on [0, oo) and m(0) = 0. 

2. Homogenization Results 
Lemma 2.1 (Barrier Functions). We define the functions ipf : M n x [0, T] — > M by 



If we chooseC > max i=li2 max (5iP ) eR n xB(0ir) \Hi(£,p)\, wherer = \\Df 1 \\ La °Q gL n ) + \\Df 2 \\L°°(R") ! 
then (<Pi,(f2 ) an d ivtiVt) are > respectively, a subsolution and a supersolution of(C e ), and 




(2.1) 



far, 0) = fa+ </?+)(•, 0) = {h,h) onR n . 



In particular, ip i < u\ < ipf on M. n x [0, T] for i = 1,2. 



Proof. We calculate that 



fai )t + Hi(~, Dtp-L ) + -fa t - <^ 2 ) 




fx{x)-h{x) 



21 



£ 
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for C > large enough as chosen above. Similar calculations give us that (fi,^) ano - 
(ipi,<pt) are 5 respectively, a subsolution and a supersolution of (C e ). By the comparison 
principle for (C e ) (see [15, 23]) we get cp7 < u\ < <pf on R n x [0, T] for i = 1, 2. □ 

Proof of Theorem 1.2. By Lemma 2.1 we can take the following half- relaxed limits 

W(x,t) := limsup%^ sup i=12 [uf](:r,t) 
w(x,t) := hminf* e ^. infj = i j2 [wf](x,t). 

We now show that W and w are, respectively, a subsolution and a supersolution of (1.3) 
in W 1 x (0, T) by employing the perturbed test function method. 

Since we can easily check W(-, 0) = w(-, 0) = / on R n due to Lemma 2.1, it is enough to 
prove that W and u> are a subsolution and a supersolution, respectively, of the equation in 
(1.3). We only prove that W is a subsolution since by symmetry we can prove that w is a 
supersolution. We take a test function <p £ C 1 (IR ri x (0,T)) such that W — <j) has a strict 
maximum at (x ,to) £ x (0,T). Let P := Dcj)(xo,t ). Choose a sequence e m — > such 
that 

W(xo,to) = limsup* maxw^ m (xo, to)- 
We define the perturbed test functions i/)^' a for i = 1, 2 and a > by 



# a 0r, !/,*) := 0(M) + £Vi(~) + 



2 



2a 2 ' 

where (1)1,1)2) is a solution of (E p ). By the usual argument in the theory of viscosity 
solutions, for every m £ N, a > 0, there exist i m>a £ {1,2} and (x mjCe , y m ,a, t m ,a) £ 
R"xK"x (0, T) such that 

El2Hn X K^ori^ m ^'^ - ^* m ' a ^' J/ '^ = U t^,S X rn,a,tm, a )-^tZ^ X ^Vm,mt m , a ) (2.2) 

and up to passing some subsequences 

("£m,cn Vm,ai ^m,a) ^ (^rai -^m, ^m) ft ^ 0, 

im,a -> i m e {1, 2} as a ->> 0, 
(z m ,t m ) ->> (x ,t ) as m ->■ oo, 
lim limw^" 1 (i m ,«,i m , a ) = W(x ,t )- 

m->-oo a— >0 ' 

Choose j m , a ,jm € {1,2} such that {i m ,a,jm,a} = {im.Jm} = {1,2}. By the definition of 
viscosity solutions, we have 

(pti^m^t tm,a)~\~Hi rn a ( , D(j)(x m ^ a , t m ,a) H ^ )H ( U i m ,a U jm,a ) ^m,a) — 0- 

(2.3) 



tvx m ,t m ) + H im (—,D(j)(x m ,t m ) + Q m ) + — {u\2 - u*2)(x m ,t m ) < (2.5) 
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Since (vi,v 2 ) is a supersolution of (E p ), we have 

H imia (—^,P + — 2 ) + [Vi ma -v jm J{——) > H(P). (2.4) 

Let a ->• in (2.3) and (2.4) to derive 

e».. 
and 

^U^, P + Qm) + (t^, - ^J(^) > H(P), (2.6) 

where Q m := lim^o (^m,a — y m , a ) / a 2 . Noting that the correctors Vi are Lipschitz continuous 
due to the coercivity of Hi, we see that \Q m \ < C for C > which is independent of m. 
Combine (2.5) with (2.6) to get 

MXm,t m ) +H(P) < H im ( — ,P + Q m ) ~ H im ( — , D((>(x m ,t m ) + Q m ) 



+ — WjZ( x m,t m ) - {4>{x m ,t m ) +e m Vj m (—))] 

~ —[ Ue C{ X m,tm) ~ {<j){Xm,tm) + £mV im { — ))} 

< cr(\P - D4>(x m ,t m )\) 

for some modulus a. Letting m — > oo, we get the result. 

We finally prove that u is Lipschitz continuous. We can easily see that / ± Mt are a 
supersolution and a subsolution of (1.3), respectively, for M > large enough. By the 
comparison principle for (1.3) we have \u(x,t) — f(x)\ < Mt for all (x,t) G M n x [0,T]. 
Moreover, the comparison principle for (1.3) also yields that 

sup \u(x, t + s) — u(x, t) | < sup \u(x, s) — f(x) | < Ms for all t, s > 0, and 

xeR n iei" 

sup \u(x + z,t) - u(x,t)\ < sup \f(x + z) — f(x)\ < r\z\ for all z G M n , t > 0. 
The proof is complete. □ 

3. Initial layers and Rate of convergence 

3.1. Inner solutions, Outer solutions, and Matched solutions. We first derive inner 
solutions, outer solutions and perform the matching asymptotic expansion method to find 
matched solutions in a heuristic mode. 

As we already obtained in Section 2, outer solutions are same as the limit u give in 
Theorem 1.2. Now we need to find a right scaling for inner solutions. We let 

<(x,t) = u\(x, -) for i = 1,2, 
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and plug into (C e ) to obtain 



[wl)t + e-Hi(p Dw{) + (w\ -w £ 2 ) = in R n x (0, T/e) 
{w £ ) t + eH 2 (^, Dw £ 2 ) + (w £ 2 -<) = in R n x (T/e), 



k ^(x, 0) = fi(x) on R n for i = 1, 2. 

We next assume that wf have the asymptotic expansions of the form 

(x, £) = Wi(x, t) + t) + e 2 Wi 2 (x, t) • ■ ■ , for z = 1, 2. 

It is then relatively straightforward to see that (wi,w 2 ) solves 

((wi) t + (w 1 -w 2 )=0 in R n x (0,oo), 
(I) i (w 2 ) t + (u> 2 - wi) = in R n x (0,oo), 

[ 0) = fi(x) on M n for z = 1, 2. 

Thus, we can compute the explicit formula for the inner solutions 

(w 1 (x,t),w 2 {x,t)) 

fi(x) + f 2 (x) fi(x)-f 2 (x) * h(x) + f 2 (x) f 2 (x)-A(x) * 
2 2 6 ' 2 2 6 

The final step is to obtain the matched solutions. We have in this particular situation 

r I *\ t /- ^ /i(^) + 
hm-u(x, r) = hm Wi(x,t) = , 

t— >0 t— >oo 2 

which shows that the common part of the inner and outer solutions is (j\ + f 2 )(x)/2. 
Hence, the matched solutions are 

/ + \ , , t\ fi{x) + f 2 (x) t fi(x) + f 2 (x) 

u(x, t) + Wi(x, -) -, u(x, t) + w 2 (x, -; 



/ + A , fl( x )~f2(x) _» , . f 2 (x)-f!(x) _2i 

*) + ^ e £ > *) + ^ e e 

= (ml(x,t),m £ 2 (x,t)), 

where mf are the functions defined by (1.4). 

As we can see, the matched solutions contain the layer parts which are essentially the 
same like the subsolutions and supersolutions that we build in Lemma 2.1. For any 
fixed t > 0, we can see that (mf (x, i), m%(x, t)) converges to (u(x, t), u(x, t)) exponen- 
tially fast. But for t = 0(e) then we do not have such convergence. In particular, we have 
(ml(x,e),m £ 2 (x,e)) converges to (u(x,t) + (/i - f 2 )(x)/(2e 2 ), u(x, t) + (f 2 - fi)(x)/(2e 2 )). 
On the other hand, the fact that ((mf)t, (m 2 )t) is not bounded also give us an intuition 
about the unboundedness of ((u\) t , (u 2 ) t )- 
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It is therefore interesting if we can study the behavior of the difference between the real 
solutions (wf,it|) and the matched solutions (mf,mf). 

3.2. Rate of convergence to matched solutions. In this subsection, we assume further 
that 

(A3) Hi are (uniformly) Lipschitz in the p- variable for i — 1,2, i.e. there exists a constant 
Cfj > such that 

\Hi(£,p) - Hi{£, q)\ < C H \p - q\ for all f £ T n and p, q £ R n . 

We now prove Theorem 1.3 by splitting M. n x [0,T] into two parts, which are M. n x 
[0,e| log e | ] and E n x [e| loge|,T]. For the part of small time W n x [0, e\ loge|], we use the 
barrier functions in Lemma 2.1 and the effective equation to obtain the results. The L°°- 
bounds of \u\ — m\\ for i = 1, 2 on R" x [e\ loge|,T] can be obtained by using techniques 
similar to those of Capuzzo-Dolcetta and Ishii [9]. 

Proposition 3.1 (Initial Layer). There exists C > such that 

-ml)(x,t)\ < Ce\\oge\ for all (x,t) £ R n x [0,e|loge|] audi = 1,2. 

Proof. We only prove the case % — 1. By symmetry we can prove the case % = 2. Let C > 
max i= i j2 niax^ ) p) 6 KT, xB (o,r) \Hi{£,p)\ be a constant, where r = ||-D/i||l°°(E") + ||-D/2||l°°(R") 
and note that u is Lipschitz continuous with a Lipschitz constant C u := (l/2)(||£)/i|| 0O + 
1 1 -C /2 1 1 oo ) • By Lemma 2.1 we have 

\ul( X ,t)-(u( X ,t) + f ^- f ^ e^)\ 

< | u ( Zj t ) - ^M±Aw | + ct = |u(x, t) - u(x, 0)\ + Ct 

< (C + C u )t<(C + C u )e\\oge\ 

for all t £ [0,e|loge|]. □ 

Proposition 3.2. Assume that (A3) holds. For T > there exists C = C(T) > such 
that 

\uf(x,t) -u(x,t)\ < Cs 1/3 for (x,t) £ R n x [e|loge|,T] and i = 1,2. 

Lemma 3.3. Assume that (A3) holds. For each 5 > and P £ M n , inere exists a unique 
solution (v(,v$) £ C^T™) 2 of 



#!(£, P + P>^(e, P)) + (1 + P) - P) = m T n , 

H 2 {t P + P>^(£, P)) + (1 + 5)^(6, P) - vHi, P) = mT». 



Moreover, 
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(i) there exists a constant C > independent of 5 such that 

P) - vt{t Q) \<C\P- Q\ for C G T", P, Q G R n and i = 1, 2; 

(ii) for each R > 0, there exists a constant C = C(R) > independent of 5 such that 

\6v?{£, P) + H(P)\ < C(R)5 for £ G T n , P G P(0, R) and i = 1, 2. 

Proof By the classical result (see [15, 23]) we can easily see that there exists a unique solu- 
tion (vf, v 5 2 ) of (E S P ) for any P G R™. We see that (v{(; P) ±C H \P - Q\/5, v 6 2 (-, P) ±C H \P - 
are a supersolution and subsolution of (El), respectively, in view of (A3). By the compar- 
ison principle for (Eg) we have 

vf(t, P) - Ch \ P ~ Q \ < vf(Z, Q) < vlii, P) + C »\ P -Q\ for e e TP, i = 1, 2, 
which completes (i). 

Let Ci(P) = max != i i2 max $6T » |P;(£,P)|. It is clear that (-C 1 (P)/6,-C 1 (P)/5) and 
(C\(P)/S, Ci(P)/5) are a subsolution and a supersolution of (Ep), respectively. Note that 

p)| < 1^,0)1 + 1^,0)-^,^)! <c(i + |P|) 

for C > max^i^^gT" \Hi(£, 0)| V Cij. Therefore, by the comparison principle again we get 

SWvfi-, P)|U~ (Tn) < d(P) < C(l + |P|). (3.1) 

Next, sum up the two equations of (E P ) to get 

P x (£, P + Dv{{£, P)) + P 2 (£, P + Dv s 2 (C, P)) < 2C(1 + |P|). 

Thus, for each P > 0, there exists a constant C = C(R) > so that 

\\Dvf(-,P) |U-(t») < C'(P) for |P| < P and z = 1,2. (3.2) 

We look back at (E p ) and take the inequalities (3.1), (3.2) into account to deduce that 

\\vi(;P)-v 6 2 (;P)\\ L o 0{Tn) < C(R) for |P| < P. (3.3) 

Let fi + := maXj =lj2 g eT n $ v i(£,yP) an d /i~ := min i=li2i £ 6 T" 8vf(£,P). Then we have 

//" < -P(P) < (3.4) 

Indeed, suppose that ji + < —H(P) and then by the comparison principle we have vf > Wt 
on T n for any solution (wi, W2, H(P)) of (E p ). This is a contradiction, since for any C 2 G R 
(wi + C 2 ,w 2 + C 2 ,P(P)) is a solution of (E p ) too. Similarly we see that yT < —H(P). 
By (3.2)-(3.4) we can get the desired conclusion of (ii). □ 
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Proof of Proposition 3.2. Let (vf(-, P), P)) be the solution of (Ef>) for P e R n . We 
consider the auxiliary functions 

$i(x, y, t, s) := u\{x, t) - u(y, s) - evf{^, - \ x ~ v\ +J f ~ s ) _ K{t + s) 

for % = 1,2, where 5 = e e and /3, 9 £ (0, 1) and il" > to be fixed later. 

For simplicity of explanation we assume that $j takes a global maximum on M. 2n x 
[e| loge|,T] 2 and let (x,y,t,s) be a point such that 

max max y, t, s) = , y, t, s). (3.5) 

i=l,2 R n x[e|loge|,T] 

For a more rigorous proof we need to add the term — 7|a;| 2 to for 7 > 0. See the proof 
of Theorem 1.1 in [9] for the detail. We first consider the case where t,s > e\ loge|. 

Claim. If < 8 < 1 - /?, then there exists M > such that (\x -y\ + \t-s\)/e p < M. 

We use y,t, s) > &x(x,x,t,t), Lemma 3.3 (i) and that u is Lipschitz continuous to 
deduce that 

< \u(x, t) - u(y, s) I + £!«?(-, *^±) - wf(* , 0) I + K\t- s\ 



2e? 



e e 



< C u (\x -y\ + \t-s\) + Csj}^-^ + K\t - s\ 
<C'(\x-y\ + \i-s\) 

for some C, C' > 0, which implies the result of Claim. 
We fix (y, s) = (y, s) and notice that the function 

(x,t) * u\{x,t) - evf(^ % -^L) - V-^^-sf _ Kt 

attains the maximum at (x,t). For a > 0, we define the function xjj by 

i;(x, e, z, t) := ul(x, t) - ev^, Z -^±) - V - ^± '* ~ ^ - |g ~ ^± |a? " ^ - Kt. 

Let ^ attain the maximum at (x a , £ a , z a , t a ) and then we may assume that (x a , £ a , z a , t a ) — > 
(x,x/e,x,t) as a — > up to passing some subsequences if necessary. By the definition of 
viscosity solutions, we have 

K H -3 h-Hi(— , s 1 1 ) + ~K -M|)(x a ,i a ) < 0, (3.6) 



and 



^ + ^-^) + (1 + <*K(£ Q , ^) - ZjL T^) > 0- (3-7) 
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Next, since ip(x a , £ Q , z a , t a ) > ip(x a ,£, a ,x a ,t a ) we get 

< e(v((U XJ1 P L ) ~ v((U ^)) < Ce^\x a - z a \ 
by Lemma 3.3 (i). Thus, \x a — z a \/a < Ce 1 ^ 9 " 13 . Send a — > and we get 

K+^ + H(^) + l(ul-ul)(x^ 

(3.8) 

Similarly we fix (x, t) = (x, t) and do a similar procedure to the above to obtain 

-K + ^ + H^^r) + C(e e + e 1 - 6 ^) > 0. (3.9) 

Combining (3.8), (3.9), and (3.5), we get 

2K <C(e 9 + e 1 ~ 9 ^). (3.10) 

Now we choose 9 = (3 = 1/3 and K = for K x large enough to get the contradiction in 

(3.10). Hence either t = —eloge or s = — 5 log 5 holds. The proof is complete immediately. 

□ 

Theorem 1.3 is a straightforward result of Propositions 3.1, 3.2. 

4. Properties of effective Hamiltonians 
4.1. Elementary properties. 

Proposition 4.1. 

(i) (Coercivity) H(P) — > +oo as \P\ — > oo. 

(ii) (Convexity) If Hi are convex in the p-variable for i = 1,2, then H is convex. 

Proof, (i) For each S > and P e R n , let (v{,v^) be a solution of (E p ) and without loss 
of generality, we may assume that i>i(£cb P) = max != i i2 ,^T" u f (£> P) f° r some £ £ T™. By 
the definition of viscosity solutions we have Hi(£ ,P) < Hi(£ ,P) + (vf — u^X^P) < 
-6vl{£o,P). We let 6 -»■ to derive that H(P) > min^i^^eT" -P)- Since i^j are 
coercive for z = 1, 2, so is iJ. 

(ii) We argue by contradiction. Suppose that H is not convex and then there would 
exist P,Q E W n such that 

2eo:=g( £±Q ) - ^); g w) >0 . (4.D 

We define the functions w { G C(T n ) so that u^f) := P) + Vi(£,Q))/2 for z = 1,2, 

where (t>i(-, P), i^Oj P)) and (t>i(-, Q), t>2(-, Q)) are solutions of (E p ) and (Eg), respectively. 
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Due to the convexity of Hi for i — 1 , 2 we have 

By (4.1) there exists a small constant 5 > such that 

#i(£, + Dwtf)) + (1 + S) Wl (C) - < - eo, 

^ tf 2 (x, + Dw 2 (0) + (1 + <5)u*(0 - < H(^^-) - so, 

and 

— ^ — + Dvi(£, — )) + I 1 + <*K(f ? —g— ) - — ^— ) > g( 2 ) - eo, 
# 2 (£, — ^— + £to 2 (f, — 2~)) + I 1 + —J - ) - —3— ) ^ g ( 2 ) - £ o- 

By the comparison principle we get 

2 < Vi(Z, 2 ) for i = 1,2. (4.2) 

Notice that (4.2) is still correct even if we replace v^, (P + Q)/2) by Vi(£, (P + Q)/2) + d 
for i — 1, 2 and for any Ci G K, which yields the contradiction. □ 

The uniqueness of the effective Hamiltonian for (E p ) and the cell problem for single 
Hamilton- Jacobi equations give the following proposition. 

Proposition 4.2. If Hi = H 2 — K, then 

H(P) = K(P) for all P eW n , 
where K is the effective Hamiltonian corresponding to K. 

Proposition 4.3. If H^ are homogeneous with degree 1 in the p-variable for i — 1, 2, then 
H is positive homogeneous with degree 1. 

Proof. Let (vi, v 2 , H(P)) be a solution of (E p ) for any P G M. n . If Hi is homogeneous 
with degree 1 in the p-variable, then (rvi,rv2,rH(P)) is a solution of (E r p) for any r > 0. 
Therefore by the uniqueness of the effective Hamiltonian we get the conclusion. □ 

Proposition 4.4. We define the Hamiltonian K as 

K(C,p) :=m a x{Hi^,p),H 2 ^,p)}. 

Let K be its corresponding effective Hamiltonian and then 

H(P) < K(P). 
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Proof. For each P e M n , there exists <p(-,P) G C°> l (T n ) such that 

K(t,P + D<p{£,P)) = K{P). 

Thus (<p(-,P),ip(;P),K(P)) is a subsolution of (E P ). We hence get K(P) > H(P) by 
Proposition 4.6. □ 

We give an example that we can calculate the effective Hamiltonian explicitly. 

Example 4.5. Let n — 1 and Hi(£,p) = \p\, H 2 (£,p) = a(£)\p\, where 

l + (| + ^)cos(2vrO >U - 
By Proposition 4.3 we have PT(P) = #(1)P for P > 0. Set 

Vl (C, 1) := sin(27rO - cos(2vrO, u 2 (£, 1) := + 77^) sin^)- 

Then we can confirm that 1), v 2 (-, !)> 1) is a solution of (Ei). Therefore P(l) = 1 

and thus, P(P) = P for P > 0. 

For any P < we have H(P) = H(-l) ■ (-P). Set 

vifa-1) ■= -(^Bin(2<) - g^cos(27rO), «a(^ 1) := "(^ + y^) sin(27r£). 

Then we can confirm that (i>i(-, —1), w 2 (-, — 1)> 1) is a solution of (E_i). Therefore H{— 1) = 
1 and thus, #(P) = -P for P < 0. Thus, we get Zf(P) = |P|. 

4.2. Representation formulas for the effective Hamiltonian. In this subsection we 
derive representation formulas for the effective Hamiltonian H(P). See [10, 20] for the 
min-max formulas for the effective Hamiltonian for single equations. 

Proposition 4.6 (Representation formula 1). We have 

H(P) = inf{c : there exists (0i,0 2 ) e C(T n ) 2 so that 

the triplet (</>i, </> 2 , c) is a subsolution o/(Ep)}. (4.3) 

Proof. Fix PeK" and we denote by c(P) the right-hand side of (4.3). By the definition 
of c(P) we can easily see that H(P) > c(P). We prove the other way around. Assume 
by contradiction that there exist a triplet (0i,0 2 ,c) G C(T n ) 2 x R which is a subsolution 
of (E p ) and c < H(P). Let i> 2 , H(P)) be a solution of (E p ) and take C > so that 
0i > Vi — C =: Ui on T n . Then since iJ* and 0j are bounded on T n , for e > small enough, 
we have 

P + Dux) + (1 + e)«i - U 2 > P + P>0i) + (1 + e)0i - <p2 
H 2 {£, P + Dv 2 ) + (1 + £)U 2 - «i > # 2 (f , P + P0 2 ) + (1 + £)0 2 - 01- 
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By the comparison principle (see [15, 23]) we deduce Ui > <pi on T n which yields the 
contradiction. □ 

If we assume the convexity on P^(£, •) for any £ £ K n , by the classical result on the 
representation formula for the effective Hamiltonian for single Hamilton- Jacobi equations 
we can easily see that 

H(P) = ^inf xn) max^C, P + D<p(£)) + t*(£, P) - v 2 & P)} (4.4) 

= ^\ Tn) p + Dm) + p) - p)] 

for any solution (v\(-, P), v 2 (-, P)) of (E p ), which is in a sense an implicit formula. For 
the weakly coupled system we have the following representation formula. 

Proposition 4.7 (Representation formula 2). If Hi are convex in the p -variable for i = 1, 2, 
then 

H(P) = inf max P + Dfctg)) + <&(0 - 4>j(€)], (4-5) 

where we take j £ {1,2} so t/iai = {1,2}. 

Lemma 4.8. Assume that Hi are convex in the p -variable. Let (v±, v 2 , H(P)) £ C(T n ) 2 be 
a sub solution o/(E p ). Setvi$(x) := /^^(rc), where p E (x) := (l/e n )p(x/£:) ancip £ C°°(lR n ) 
be a standard mollification kernel, i.e., p > 0, supp p C P(0, 1), anc? L„ p(x) dx = 1. Then, 
(vis, v 2 s, H(P) + U)(5)) is a subsolution of (E p ) for some modulus u. 

Proof. Note that in view of the coercivity of H iy V{ are Lipschitz continuous and (v±, v 2 , H(P)) 
solves (E p ) almost everywhere. Fix any £ £ T™. We calculate that 

H(P) > p s *{H 1 (;Dv 1 (-)) + (v 1 -v 2 ))(C) 

= / Ps(£-v)( H i(v, D Mv)) + ( v i- v 2)(v))dv 

> I Ps(C - v)(Hi(£, DvM) ~ drj + (v 15 - u M )(0 

> P s * D Vl {£)) + (v 15 - v 2S m - u(S) 
= Dv ls (0) + (vu - v 2S )(0 - u(5), 

where the third inequality follows by using Jensen's inequality for some modulus u. □ 

Proof of Proposition 4.7. Let c(P) be the constant on the right-hand side of (4.5). Noting 
that for any ((f>i,<f) 2 ) £ C^T") 2 

H^, P + D(f>i(£)) + - < max [H^, P + + - &)(0] == ^ 

for every £ £ T n . By Proposition 4.6 we see that H(P) < a M2 for all (0i, 2 ) £ C^T™) 2 . 
Therefore we get H(P) < c(P). 



18 H. MITAKE AND H. V. TRAN 

Conversely, we observe that by Proposition 1.1 (v±s(-, P), V2s(-, P), H(P) + cu(5)) G 
C 1 (T n ) 2 x M is a subsolution of (E p ). Therefore, by the definition of c(P) we see that 
c(P) < H(P) + u(5). Sending 5 — > yields the conclusion. □ 

If Hi are convex in the p- variable, then there is a variational formula of solutions of the 
initial value problem and the cell problem as stated in Introduction. Therefore, naturally 
we have the following variational formula 

- + 0O 



/-t-oo 
e~ Ss (-P ■ f](s) + L„ {s) ( V (s), - v (s))) ds 

(-P ■ fj{s) + L u{s) (r](s), -t?(s))) ds 



1 



= — lim - inf E, 

t— >oo t V LJq 

where the infimum is taken over r\ G AC ([0, +oo), IR n ) such that r/(0) = x and Ej denotes 
the expectation of a process with u(0) = i given by (1.2). 

Remark 4.9. When we consider the nonconvex Hamilton- Jacobi equations, in general we 
cannot expect the formula (4.5). Take the Hamiltonian 

Hifop) :=(|p| 2 -l) 2 fori = 1,2 (4.6) 

for instance. In this example if we calculate the right-hand side of (4.5) with P = 0, then 
it is 0. But we can easily check that H(0) = 1, since in this case we have solutions (which 
are constants). 

The following formula is a revised min-max formula for the effective Hamiltonian for 
nonconvex Hamilton- Jacobi equations. 

Proposition 4.10. We have 

H(P)= inf max sup [#*(£, P + p) + (0, - ^)(0], ( 4 -7) 

where if D + <pi(£) = 0, then we set sup peD+(# ,. ( ^ [H^, P + p) + - 0j)(O] = _ °° 
convention. 

We notice that if Hi are given by (4.6), then the right-hand side of (4.7) with P = is 

I. 

Proof. The proof is already in the proof of Proposition 4.7. We just need to be careful for 
the definition of viscosity subsolutions. Indeed, let c be the right-hand side of (4.7) and 
noting that for any (fa, fa) G C°> l {T n ) 2 , £ G T n , and q G D + fa(£), 

Hi{£, P + q) + (fa- fa)(0 < max sup [#*(£, P + p) - (^ - ^)(01 = : ^i,^- 

C6T» l t=l,2 p6D+ ^.( fl 

Thus, H(P) < a^ 2 for all (fa, fa) G C '^!™) 2 by Proposition 4.6. Therefore, H(P) < c. 
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1!) 



Conversely, there exists a viscosity subsolution (t>i(-, P), v 2 (-, P), H(P)) G C 0,1 (T n ) 2 x 
of (E p ). By the definition of viscosity subsolutions we have 

Hi{t P + p) + (vi- Vj)(£) < H{P) for all (GT n and p G D+v^). 

Thus, 

max sup [Hi(£, P + p) + {v, - Vj )(£)] < H(P), 



which implies c < H{P). 



□ 



4.3. Flat parts and Non-flat parts near the origin. In this subsection, we study the 
results concerning flat parts and non-flat parts of the effective Hamiltonian H near the 
origin. We first point out that there are some cases in which we can obtain similar results 
to those of Concordel's results for single equations. We present different techniques to 
obtain these results , namely the min-max formulas, and the construction of subsolutions. 
In this subsection, we only deal with the Hamiltonians of the form Hi(£,p) = \p\ 2 — V$(£), 
where Vj G C(T n ) for i — 1,2 unless otherwise stated. 

Theorem 4.11. Assume that V { > in T 1 and {V { = 0} =: U { C T n for % = 1,2. We 

assume further that U\ D U 2 ^ and there exist open sets W\,W 2 in T n , and a vector 
qeR n such that U{q + W 2 ) <s (0, l) n and 

[/iU(7 2 c^C W 2 and dist(Wi, dW 2 ), dist(E/i U U 2 , dW x ) > 0, (4.8) 

then there exists 7 > such that H(P) = for \P\ < 7. 





Fig. 4.1. The figure of Ui, W{. Fig. 4.2. The graph of p in case n = 1. 

Proof. Without loss of generality, we may assume that q = 0. Take £0 G U\ H L^. By 
Proposition 4.7 we have 

H{P) > inf max[|P + J D^o)| 2 -^(eo) + ^(eo,P)-^(eo,P)] 

> inf max [</>i(£ , P) - </?i(£o, P)] = 0. 
(vi,va)GC 1 (T») 1=1,2 Lr ' rjv 
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Now, let d := min{dist(Wi, dW 2 ), dist(E/i U U 2 ,dW 1 )} > 0. There exists e > such 
that 

Vi(0 > e for ieT"\^, i = 1, 2. (4.9) 
We define a smooth function ip on T n such that 



(p(g) = -P . £ on Wi, <^(0 = on T n \ W 2 
1 



\D<p\ < on T n . (4.10) 



Notice that 



= -^(0 < 0, on W u 

< - £o, on T™ \ Wt. 

Thus, |P + Dip{£) | 2 - < on T n provided that |P| < d^/C =: 7. We hence have 
that (ip,(p,0) is a subsolution of (Ep) for |P| < 7. Therefore H(P) < for |P| < 7 by 
Proposition 4.6. □ 

Remark 4.12. (i) In fact, the result of Theorem 4.11 still holds for more general Hamilto- 
nians 

where P G C(T n x R n ) and G C(T n ) are assumed to satisfy 

(a) the functions p h-» Pj(£,p) are convex and Pj(£,p) > Pj(£,0) = for all G 
T n x R n , 

(b) Vi, V2 satisfy the conditions of Theorem 4.11. 

(ii) Notice that the assumptions V, > and {Vi = 0} 7^ for z = 1, 2 are just for simplicity. 
In general, we can normalize V, by Vi — ming 6 T« to get back to such situation. 

(iii) From the proof of Theorem 4.11 we have 

P) = u 2 (£, P) for all £ G {V x = 0} n {V 2 = 0} 
for any solution (vx(- , P) , v 2 (- , P)) of (Ep). 

(iv) By Proposition 4.4 we can give another proof to Theorem 4.11 as follows. In this case, 
we explicitly have 

K(£,p) = max{H 2 - V(C), \p\ 2 - V 2 (C)} = \p\ 2 ~ V(C) 

where V(£) = min{Vi(£), ^(0>- Note that V > and {1/ = 0} = {Vi = 0} U {V 2 = 0}. 
Hence, we can either repeat the above proof for single equations to show that H(P) < 
K(P) = for |P| < 7 or we can use Concordel's result directly. 
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Notice that condition (4.8) is crucial and plays an important role in the construction of 
the subsolution (ip, ip, 0) of (Ep) and could not be removed in the proof of Theorem 4.11. 
We point out in the next Theorem that there are cases when (4.8) does not hold, then the 
flatness near the origin of H does not appear. 

Theorem 4.13. Assume that Vi > in T n and 

{Vi = 0} = {V 2 = 0} = {£ = (6, • • • G TT : & = 1/2 for j > 2} =: K. (4.11) 
The followings hold. 

(i) There exists 7 > such that H(P) = \P\\ 2 provided that \P'\ < 7 for any P = 
(Pi,P') 6Rx R n -\ 

(ii) H(P) > I Pi 1 2 for all P e R n . 

Proof. Firstly, we prove that H(P) < \P\\ 2 provided that \P'\ < 7 for some 7 > small 
enough by using exactly the same idea in the proof of Theorem 4.11. We build a function 
f(0 — <f(£,2, ' ' ' ;£n)G C 1 (T n ), which does not depend on £ 1; so that 

n 

X;|Pi + ^(0l a - W<0on T n 
i=2 

for i — 1, 2 and for |P'| < 7 with 7 > small enough. Thus 

\P + D<p(Z)\ 2 -Vi(£)< |Pi| 2 on T n 

for i = 1, 2. By Proposition 4.6 2F(P) < |Pi| 2 . 

We now prove that H(P) > \Pi\ 2 . For each £0 G -fT, we have in view of (4.4) 

H(P)= inf max[\P + D<p(t)\ 2 -V 1 (t)+v 1 &P)-v 2 (t,P)} 
> inf [\P + D< f (^)\ 2 + v 1 (Co,P)-v 2 (Co,P)], 

and similarly 

H(P) > inf ' [\P + D^ )\ 2 + v 2 (£ , P) - Vl (Zo, P)]. 
Take an arbitrary function ip G C 1 (T n ) and observe that 

/ |p + ^(OI 2 ^i 

Jk 

> [ |p 1+ ^(oi 2 ^i= / \Pi\ 2 +\^m 2 + 2^(0^1 

> / |Px| 2 + 2^(0^1 = |Pi| 2 . 

Thus, it is clear to see that H{P) > \Pi\ 2 , which implies the result. □ 
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The above two Theorems describe several examples that we can obtain similar results of 
the flat part or non-flat part of H to those of single Hamilton-Jacobi equations in [11, 12]. 
Indeed, the structures on the potentials Vi for % — 1, 2 are very related in such a way that 
we obtain the shape of H like for single equations. We rely on the idea of building the 
subsolutions (<p, ip, H(P)) of (E p ) where (p = ip, which does not work in general cases. 

Next, we start investigating the properties of H in some cases where the structures of 
the potentials Vi for i — 1, 2 are widely different and in general we cannot expect H to have 
simple properties. The next question is that: Can we read of information of the effective 
Hamiltonian in the case where {V\ = min^ gT n Vi(£)} D {V2 = min^gT™ V 2 (£,)} — 0? 

Theorem 4.14. Letn = 1 and assume that for Eq > small enough the following properties 
hold. 

4 12 3 13 

(a) {Vi = 0} = [-, -], {Vi = -e } = [0, 1] \ (-, -), and -e < V x < on T for 

lb lb lb lb 

some Eq > 0. 

7 9 b 10 

(b) {V 2 = 0} = [-, -], {V 2 = 2} = [0, 1] \ (-, -), and < V 2 < 2 on T. 

There exists 7 > such that H(P) = for \P\ < 7. 
Lemma 4.15. We have 

H(P) > -lmm(V 1 + V 2 )(0- 
2 ?gT" 

Proof. Sum up the two equations in (Ep) to get 

\P + Dv ± \ 2 + \P + Dv 2 \ 2 -Vi-V 2 = 2H(P), 

which implies 2H(P) > — (Vi + V 2 )(l;)) for a.e. ^ G T n , and the proof is complete. □ 

Proof of Theorem 4.14. Noting that min eeT (Vi+V 2 )(0 = 0, and {Vi = -e}n{V 2 = 0} = 0, 
we have H(P) > by Lemma 4.15. We construct a subsolution (<£>, -0,0) of (E p ) for \P\ 
small enough. Let 

v lb lb ; v lb lb ; v lb lb ; 

Let P < for simplicity. We define the functions </?, ip by 

for x e W2 
for f e T \ 1V 3 

and \D(p\ < C X \P\ for some d > 0, < ip < -P ■ £ on [0, 1] and 

-P • f for £ G Wi 

C 2 £eT\PV 2 

for some C 2 G (1/64, 1), |P + < 1, and ifj > -P ■ £ on [0, 1]. 



V(0 := 




^(0 



HOMOGENIZATION OF WEAKLY COUPLED SYSTEMS OF HJ EQUATIONS 23 

We have 

\p + d^)\ 2 - K(0 + v(0-^(0 
< (V(0-^(0<o if£ew 2 

~ \2{Cl + 1)|P| 2 + e + \P\ -C 2 iii&T\W 2 . 
If \P\ and e are small enough, then \P + | 2 — Vi(£) + </?(£) — VKO < on T. Besides, 

|P + J D^)| 2 -F 2 (£) + ^(£) -</?(£) 

~jl-2 + C7 2 -0<0 if£eT\Wi. 




W 3 W 2 
Fig. 4.3. The graph of Fig. 4.4. The graph of ip. 

Thus i/), 0) is a subsolution of (Ep), and the proof is complete. □ 



Remark 4.16. It is worth to notice that 

H(0)^~mm(V 1 + V 2 )(x) 

in general. Indeed, set 

Vi(0 = 4tt 2 sin 2 (27rO + cos(2<) - sin(2<), 
y 2 (£) = 4tt 2 cos 2 (2<) + sin(27r£) - cos(2tt£). 

If we check that (cos(27r£), sin(27r£), 0) is a solution of (E ), then we realize H(0) = 0. In 
this case 

if (0) = + -2tt 2 = --(V 1 + y 2 )(0 for all £ G T. 

In Theorem 4.11 the fact that U(R n \ (E/i U C/ 2 )) is connected, where Ui = {Vi = 0} plays 
an important role in the construction of subsolutions as stated just before Theorem 4.13. 
In the next couple of Theorems we make new observations that we can get the flat parts 
of effective Hamiltonians even though II(M n \ (U\ U U 2 )) is not connected. 
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Theorem 4.17. Let n = 1 and assume V\ = 0, V 2 > on [0, 1] and {V 2 = 0} = {1/2}. 
Then there exists 7 > such that H(P) = for |P| < 7. 

Sketch of Proof. The proof is almost the same as the proof of Theorem 4.14 but let us 
present it here for the sake of clarity. Since min^in^i + V2XO = 0, we have H(P) > 
by Lemma 4.15. Let 

There exists M G (0, 1) so that 

V 2 (0 > M for i£ W v 

Assume P < for simplicity. We now construct the functions <p, i/j so that -0,0) is a 
subsolution of (E p ) for small \P\, which implies the conclusion. Take \P\ < M/4 first. We 
define the functions <p, ip by 

-P ■ £ for x G W 2 

for £ G T \ W 3 

and |P></?| < C X \P\ for some C x > 0, < y? < -P ■ £ on [0, 1] and 

\c 2 £gT\W 2 

for some C 2 G (M/128, M/2), |P + Dijj\ < M/2, and ^ > -P • £ on [0, 1]. 
We have 

IP + W + rtfl-rtf)^ - * * ,f(a ~ W > 

mi\ t \ 2 (Cf + l)\P\ 2 + \P\-C 2 if£eT\W 2 . 

If |P| is small enough, then |P + Dif(C,)\ 2 + — V>(£) < on T. Besides, 

|p + ^(OI 2 -v 2 (0 + ^(0-^(0 

if £ e Wi 

M 2 M 
M + C 2 -0<— -M+— <0 if £ G T \ W\ . 

Thus (ip,ip,0) is a subsolution of (Ep), and the proof is complete. □ 

We can actually generalize Theorem 4.17 as following. 

Theorem 4.18. Assume that Vi = 0, V 2 > and there exist an open set W in T n and a 
vector q el" such that U(q + W) <s (0, l) n and ^ {V 2 = 0} C W. Then there exists 
7 > such that H(P) = for \P\ < 7. 

The proof of this Theorem is basically the same as the proof of Theorem 4.17, hence 
omitted. The following Corollary is a direct consequence of Theorem 4.18 
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Corollary 4.19. Assume that Vi, V 2 > and t/iere exzst an open set W in T n and a vector 
qeR n such that U(q + W) m (0, l) n and 

{v x = o} n {V 2 = 0} c {V 2 = 0} c w. 

Then there exists 7 > such that H(P) = for |P| < 7. 

The result of Corollary 4.19 is pretty surprising in the sense that flat part around of 
H occurs even though we do not know much information about V\. More precisely, we 
only need to control well {V 2 = 0} and do not need to care about {Vi = 0} except that 

{Vi = 0} n {v 2 = 0} ^ 0. 

Finally, we consider a situation in which the requirements of Theorem 4.18 and Corollary 
4.19 fail. 

Theorem 4.20. We take two potentials V % : T — > [0, 00) such that V 1 are continuous and 
{V 1 = 0} = {yoi} for some yoi € T for i — 1,2. Assume that VI (£1,62) = ^(^i) an ^ 
^2(61,6) = V 2 {£ 2 ) for (£i,£ 2 ) G T 2 . r/ien tnere exists 7 > sncn that H{P) = /or 
\P\ < 7- 





{^2 = 0} 




L=0} 



yoi 

Fig. 4.5. The figures of {Vj = 0} 

Proof. By using Theorem 4.17, for P = (P l7 P 2 ) with \P\ small enough, there exist two 
pairs {(f u ipi) G C°^(T) 2 for i = 1,2 such that 

Pi + v?'i(£i)i 2 - ^(6) + ^(Ci) - Mti) = 0, 
A+^i(ei)i 2 +^i(ei)-^(6) = o 

and 

P2 + <p 2 (&)\ 2 + M&)-M&) = u, 
p 2 + ^ 2 (6)| 2 - v\&) + ^(6) - ^2(6) = 

Now let = <M£i) + ^2(6), ^(£i,£ 2 ) = ^i(a) + ^(£2) for (6,6) G T 2 . For 

P = (Pi,P 2 ) with |P| < 7, we easily get that (t>i,t> 2 ,0) is a solution of (Ep), which means 

pT(p) = 0. □ 
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5. Generalization 
In this section we consider weakly coupled systems of m-equations for m > 2 



j m 

+ Hi(-, Duf) + - V Cji(«? - O = in R™ x (0, T) for i = 1, . . . , m, 

P P ' ^ 



with 



0) = on R n for i = 1, . . . , m, 



where Cy are given nonnegative constants which are assumed to satisfy 



£ 

i=i 



Cjj = 1 for all i = 1, 



m. 



(5.1) 



Set 





' en 


• • • Ci m \ 




f «! 








\ 


K := 












, and f := 








\ c m i 


Cmm y 










^ /m 


/ 



Then the problem can be written as 

/ Hx(x/£,Dul) \ 



+ -(I-K)u £ = inR n x(0,T), 



(5.2) 



y H m (x/£,Du e m ) J 
Lu e (-,0) = f onR n , 
where / is the identity matrix of size m. We obtain the following result. 

Theorem 5.1. The functions u\ converge locally uniformly to the same limit u in R n x 
(0, T) as e —> for i — 1, . . . , m and u solves 

u t + H(Du) = in R n x (0, T) 
u(x,0)=J(x) onW 1 , 

where H is the associated effective Hamiltonian and 



1 m 

J(x) : 

i=i 



We only present barrier functions which are generalizations of the barrier function in 
case m = 2 defined by (2.1) in Lemma 2.1. Set 

w ± (x,t) := (7±Ct)j + g £ (x,t), 

where C is a positive constant which will be fixed later, j := (1, . . . , 1) T and 

g £ (x,t) := [e^ K -Vh\(x), h(x) := f (x) - J(x)j. 
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Since we assume (5.1), we can easily check that the Frobenius root of K, i.e., the maxi- 
mum of the eigenvalues of K, is 1 and moreover j is an associated eigenvector. Moreover 
by the Perron-Frobenius theorem we have 

Lemma 5.2. There exists 5 > such that |e^ /<_/ ^h| < e _<s '|h| provided that h • j = 0. 

See [24, Lemma 5.2] for a more general result. 

Proposition 5.3. The functions w ± are a subsolution and a super solution of (5.2) with 
w ± (-,0) = f on M. n , respectively, if C > is large enough. 

Proof. It is easy to check w ± (-, 0) = f on IR n . Note that 

= -(K - I)g £ and |Dg| < Ce~% . 
dt e ( 1 

Thus, we can check easily that are a subsolution and a supersolution of (5.2), respec- 
tively, if C > is large enough. □ 

By a rather standard argument by using the perturbed test functions we can get Theorem 
5.1 as in the proof of Theorem 1.2. 

6. Dirichlet Problems 

In this section we consider the asymptotic behavior, as e tends to 0, of the viscosity 
solutions (uf , w|) of Dirichlet boundary problems for weakly coupled systems of Hamilton- 
Jacobi equations 



u\ + #x(p Du\) + -(«?- 4) = in Q, 
u\ + H 2 (-,Dv,2) + -{u% - u\) = in Q, 



^ u\(x) = gi(x) on dQ, 

where Q is a bounded domain of M. n with the Lipschitz boundary, the Hamiltonians Hi G 
C(R n x W l ) are assumed to satisfy (Al)-(A2) and gi G C(dQ) are given functions. 

Concerning the Dirichlet problem, classical works required continuous solutions up to 
the boundary to satisfy the prescribed data on the entire boundary. This can be achieved 
for special classes of equations by imposing compatibility conditions on the boundary data 
or by assuming the existence of appropriate super and subsolutions. However, in general, 
we do not expect that there exists a (viscosity) solution satisfying the boundary condition 
in the classical sense. After Soner studied the state constraints problems in terms of PDE, 
the viscosity formulation of Dirichlet conditions was introduced by Barles and Perthame [3] 
and Ishii [22] . In this paper we deal with solutions satisfying Dirichlet boundary conditions 
in the sense of viscosity solutions. 
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Theorem 6.1. Let be the solution o/(D £ ). Then uf converge locally uniformly to 

the same limit u on Q as e — > for i — 1,2 and u solves 

u + H(Du) = inn, 

V _ (6.1) 
u(x, t) = g on <9f2, 

where ~g := min{gi, g 2 } on dVt. 

Lemma 6.2. If(u\,u £ 2 ) is a bounded uppers emicontinuous subsolution of (D £ ) , thenuf(x) < 
gi(x) for all x G dVt and i = 1,2. 

Proof. Fixxo G dQ. Choose a sequence {xfcjfceN C M n \f2 such that \xq — Xk\ = 1/k 2 . Define 
the functions </>i : f2 — > IR by (f>i(x) := u\ (x) — k\x — Xk\- Let r > and G B(xo, r) fl f2 
be a maximum point of 0i on B(x ,r) D f2. Since 0i(^) > 0i(x o ), we have — x k \ < 
M i(£fc) — u l( x o) + ^l^o — ^fcl < Cj where C > is a constant independent of k. Thus, 
io as I; ^ 00. Moreover, noting that u\(xq) < liminffc^ 00 (M^ + &|xo — Xk\) < 
limsup^^M^^fc) + limsupfc^ k\x - x k \ < u\{xq), we get tif(Cfe) -> wf(x ) as -> 00. 
By the viscosity property of u\ , we have 

+ #i(^,p fc ) + -(«?(&) - «2(&)) < or (6.2) 

where = &(£fc — ^fcVl^fc — #fc|- Noting that \pk\ = ^, by (Al), we see that the left-hand 
side of (6.2) is positive for a sufficiently large k G N and then we must have u\{^k) < 9i{£,k)- 
Sending k — >■ 00, we get wf (x ) < gi^o)- Similaly, we get m|(x ) < #2(^0) 011 <9fi. D 

Lemma 6.3. The families {u e } £> o are equi-Lipschitz continuous in Q for i = 1,2. 

Proo/. Set M := max 4=li2 (||^(-, 0)|| L =c (Rn) + ||^|U-(^))- Then (— M, — M) and (M, M) 
are a subsolution and a supersolution of (D £ ), respectively. By the comparison principle 
for (D £ ) we have \uf\ < M. Adding two equations in (DJ we get 

u\ + u £ 2 + #i(~ Du\) + H 2 (^, Du £ 2 ) = 

for almost every x G f2, which implies that \Duf\ < M in the sense of viscosity solutions 
for some M' > 0, which is independent of e. □ 

Proof of Theorem 1.2. By Lemma 6.3 we can extract a subsequence {ej} converging to 
so that Ui converges locally uniformly to Ui G C(Q) for i = 1,2. By usual observations, 
we get that u± = u 2 =: u. Since (6.1) has a unique solution, it is enough for us to prove 
that u is a solution of (6.1). 

We only prove that u is a supersolution of (6.1), since in view of Lemma 6.2 we can 
easily see that u is a subsolution of (6.1). 
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Let G C 1 (r2) be a test function such that u — <fi takes a strict minimum at xq G 
Q. We only consider the case where xq G dQ, since we can prove by a similar way to 
the proof of Theorem 1.2 in the case where Xq G f2. It is enough for us to prove that 
u(x ) + H{D<P(xq)) > provided that (u — g)(x ) < 0. 

Let (vi,v 2 ) be a solution of (Ep) with P := D(f>(x ). We consider 

m e := min min(uf(x) — — evA— )). 

i€{l,2} xG fi £ 

Pick 2 e G {1,2} and x e G f2 so that m e = u £ e {x £ ) — <fi(x £ ) — evis(x £ /e). Also choose 
j £ G {1,2} such that {i £ ,j £ } = {1,2}. We only consider the case where x e G dQ again. 
Since Ui £ converges to u locally uniformly on Q, ev{-/e) converges to uniformly on Q as 
e —7- and u — <p takes a strict maximum at Xq, we see that x £ — > Xq as e — > 0. Thus, if e 
is small enough, then we may assume that {u% — g^i^x 6 ) < 0. 
For a > we define the function <3? Q : f2 x M. n -> E by 



?/ 1 1 

$ a (x, y) := - 4>(x) - ev^ (-) + —\x - y\ 2 + -|x - £' 



e |2 



2a 2 ' 1 2' 

Let $ a achieve its minimum over Q x R n at some (x^, y^). Since we may assume by taking 
a subsequence if necessary that x £ a — > x £ as a — > 0, we have 

(uie — gis)(x £ a ) < for small a > 0. 

Therefore, by the definition of viscosity solutions, we have 

u ie + H t ^,D^(x £ a )-p £ a - (x £ a - x £ )) + - e {u% - u £ s )(x £ a ) > 0, 

where p £ a := (x £ a — y £ a )/a 2 . Also, we have 



fr,(£ p - p £ J + - ^)(f ) < ^(P), 

since (^1,^2) is a solution of (Ep). 

A priori Lipschitz estimate implies \p e a \ < C for some C > which is independent of a 
and £. Without loss of generality, we may assume that p £ a — > p £ by taking a subsequence 
{«j} converging to if necessary. Send a — > in the above inequalities to obtain 

u%(x £ ) + H ie ( — ,D(j)(x £ ) -p £ ) + -{u%(x £ ) - u £ (x £ )) > 0, 
e e 

P-p £ )+ - ^e(-) < tf(P). 

Noting that u%(x £ ) — <p(x £ ) — ev^^x 6 /e) < u E e (x e ) — <p(x £ ) — evje(x £ /e), we get that 

uUx £ ) + H(P)>H i ,(j,D(l)(x £ )-p £ )-H i ,(^-,P-p £ ) > - a (\D<P(x £ )-P\). 
Sending e — > yields the conclusion. □ 
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In order to explain the relation between (D e ) and the exit-time problem in the optimal 
control theory, we assume that the Hamiltonians Hi are convex in the p- variable henceforth. 
We next define the associated value functions, which give us an intuition about the effective 
boundary datum ~g in Theorem 6.1. 

For e>0we define the functions itf : Q — > M by 

u\{x) := inf {Z, ( £ e-L^pQ, -f)(s)) ds + e~ T g^(r)(v(r)))}, (6.3) 

where the infimum is taken over n G AC ([0, oo), Q) such that 77(0) = x and r G [0, 00] such 
that 77 (r) G dQ and if r = 00, then we set e~°° := 0. Here denotes the expectation of a 
process with v £ {f)) = i, where v e is a {1, 2}-valued continuous-time Markov chain given by 
(1.2). 

Theorem 6.4. Assume that the functions u\ given by (6.3) are continuous on Q. Then 
the pair (wf,w|) is a solution o/(D e ). 

The proof of Theorem 6.4 is given in Appendix. See [3, 22] for single equations. The 
value functions defined by (6.3) give us an intuitive explanation of the reason why the 
boundary datum ~g of the limit solution u is the minimum of for % = 1, 2. If we send e to 
0, then the switching rate becomes very fast but it does not really affect the exit time as we 
can choose to stay in Q as long as we like. And hence, we can control the exit state in such 
a way that the exit cost is the minimum of two given exit costs On the other hand, when 
we consider the value function (1.1) associated with the initial value problem, we cannot 
control the terminal state and also the timing of jumps, which are only determined by a 
probabilistic way given by (1.2). This is the main difference between Dirichlet problems 
and initial value problems and the reason why the effective Dirichlet boundary value and 
the effective initial value are different. 

7. Appendix 

We now prove Theorems 7.1, and 6.4 by basically using the dynamic programming 
principles, which are pretty standard in the theory of viscosity solutions. Throughout 
this section we always assume in addition to (Al), (A2) that p 1— y Hi(£,p) are convex for 
i = 1,2. 

Theorem 7.1 (Verification Theorem). Assume that the functions u\ given by (1.1) are 
continuous on IR n x [0,T]. Then the pair (ttf,w|) is a solution o/(C e ). 

Let e — 1 for simplicity in what follows. By abuse of notations we write (^1,^2) for 
KX)- 
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Proposition 7.2 (Dynamic Programming Principle). For any x G IR n , < h < t and 
2 = 1,2 we have 

Ui(x,t) = inf { E i(y L u ( s )(r](s),-fi(s))ds + u v ( h )(r](h),t - hj^Y (7.1) 
where the infimum is taken over rj G AC ([0, /i],M n ) with r](0) = x. 

Proof. We denote by Vi(x,t;h) the right-hand side of (7.1). Let rj be a trajectory in 
AC([0,t],R n ) with 77(0) = x and v be a process with z/(0) = 2 which satisfies (1.2). Set 
fj(s) := r](s + /i) and i>(s) := v(s + /i) for s & [Q,t — h]. We have 





/ £„( 8 )(t7(s), 
'0 


-r}(s))ds + . 


k(t)fo(*))) 




/ ^(,)(^(s), 
/o 


—f](s)) ds + 


J L v{s) {r){s), -rj(s)) ds + U( t )(v(t))) 


= E<( 


/ ^)(^(s), 


-r)(s))ds) - 


|-E^y L Hs) (fj(s),-fj(s)) ds + f i )(t-h)(fj( t 


>Et( 


Jo 


—f](s)) ds + 




> fj(x. 


,t;h), 







which implies Ui(x,t) > Vi(x,t; h). 

Let 5i G AC([0,/i],M n ) and 5 2 G AC([0,t - /i],M n ) be trajectories with S^h) = S 2 (0) 
and 5i(0) = x. Set 

^ : ~ I S 2 (s-h) for all s G [M]. 
Let z/ be a process with z/(0) = 2 which satisfies (1.2). Note that 



h ft—h 

Lv( s ){5i{s), -6i(s)) ds + / L u{s+h) (5 2 {s), -5 2 {s)) ds + f u {t)(S 2 (t - h)) 

Jo 

t 

Lv(s){v(s), -fj{s)) ds + fu(t){v{t))- 



We have 



rh ft—h 

E *(y L u{s - ) (5i(s),-6i(s))ds + J L u{s+h) (5 2 {s),-5 2 {s))ds + f u{t) {5 2 {t - h)) 
= Ei {f o L ^s)(v(s), -V(s)) ds + f v( t) (r/(t)) 

>Ui(x,t). 
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Take the infimum on all admissible 62 to obtain 

E i(y L^ a )(Si(s) i -Si(8))d8 + u v ( ih )(S 1 {h) i t-h)j >Ui(x,t), 

which implies Vi(x, t; h) > Ui(x,t). □ 

Proof of Theorem 7 .1 . It is obvious to see that (1x1,1x2) (-,0) = (/lj/b) on R™. We first 
prove that u\ is a subsolution of (Ci). We choose a function £ C 1 (lR n x (0, T)) such that 
ui — has a strict maximum at (xo,t ) £ R n x (0, T) and (tti — <f)(xo,to) = 0. 
Let /i > 0. By Proposition 7.1 we have 

ui(x ,t ) < E i(^ L u{s )(r](s),-i)(s))ds + u l , w (r](h),t - hfj (7.2) 

for any 77 £ AC ([0, R n ) with 77(0) = xq £ R n and r)(0) = q £ R n . Since 1/ is a continuous- 
time Markov chain which satisfies (1.2), the probability that u(h) = 2 is C\h + o(h) and 
the probability that u(h) = 1 is 1 — (c\h + o(h)). By (7.2) we obtain 

<f>(%o,to) = ui(x ,t ) 

< (1 - Ci/i - o(/i)) ( y L x (^, -7)) ds + Ut(T](h),t - /i)) 

+ (ci/i + o(/i)) ( J L 2 (r), -rj) ds + u 2 (r](h),t - /i)) 

< / L 1 (r],-ri)ds + (j)(r](h),tQ-h) 
Jo 

+ (ci/i + o(/i)) ( / L 2 (t7, -77) (is + u 2 (rj(h),t - h) - Li(r), -rj) ds - ui(r)(h),t - hfj . 



Thus, 

(j>(r}(0),to)-(l>{ri(h),to-h) 
h 

< ~ y Li(t7, -77) ds + (ci + ^)(u 2 (v(h),t - h) - m^h),^ - h)) 

+ (ci + ^)(y L 2 (r),-r))ds- y L r (ri,-rj)ds). 
Sending /i — >■ 0, we obtain 

<Pt( x o,t ) + D(p(x ,t ) ■ (-q) < L ± (x , -q) + ci(u 2 - «i)(a?o,*o) for a11 Q e R n , 
which implies (j)t(x ,t ) + H(x , D(j)(x ,t )) + c±(u 2 - Mi)(x ,to) < 0. 
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Next we prove that u 1 is a supersolution of (Ci). We choose a function <fi G C 1 (IR n x 
(0, T)) such that u±— (f> has a strict minimum at (xo, to) G R n x (0, T) and (ui—(f))(xo, to) = 0. 
Take h, 8 > 0. By Proposition 7.1 we have 

Mi(xo,*o) + S >Ei {J Q L ^(s)(V5(s), -risis)) ds + u u{h )(r] 5 (h),to - &)) (7.3) 

for some r]s G AC ([0, /i],M n ) with 77,5(0) = Xo- Since v is a continuous-time Markov chain 
which satisfies (1.2), by a similar calculation to the above we obtain 

<t>(x , t ) + 5 = ui (x , t ) + 5 
> / Lx(r),-ri)ds + <l>(7)(h),tQ-h) 



+ {c 1 h + o(h))[ I L 2 (r),-f))ds + u 2 (r)(h),to- h) - I L^r), -77) ds - «i (77(h), t - h) 



Thus, 



5 1 f h d(f)(r]s(s),to- s) 

h > hj ds + L ^s,-Vs)ds 

+ (ci + ^-r-)(u 2 (7]s(ti),t - h) - ux(r) S (h),t - h)) 



o{h) 

+ (ci + 



1 



(f>t(r]s(s),t - s) - D(f>- (-rjs(s)) + £1(775, -r/'a) 



+ (ci + ^)(u 2 (vs(h),t -h)- Ul (r]s(h),to - h)) + O(h) 
1 f h 

~hj -( < t>t(V6(s),t -s) + H 1 {r ]s {s),D<f)))ds 

+ (ci + ^-)(u 2 (r]s(h),t -h)- ui(rj S (h),t - h)) + 0(/Y). 

We finally set 5 = h 2 and let /i — >■ to yield the conclusion. □ 
By a similar argument to the proof of Proposition 7.2 we can prove 

Proposition 7.3 (Dynamic Programming Principle for (6.3)). For any x G M n , h>0 and 

i = 1,2 we have 

Ui(x) =inf A e- s L u{s) (^,-fi(s))ds 

+ l {h<T} e~ h u u{h) {rj(h)) + l {h > T] e~ T g u{T) {rj(T)) ) J, (7.4) 
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where v with z/(0) = i is a {1, 2}-valued continuous-time Markov chain which satisfies (1.2) 
and the infimum is taken over 77 G AC ([0, h], Q) such that 77(0) = x and r G [0, h] such 
that r](r) G dft. 

Proof of Theorem 6.4. We only prove in what follows that it, satisfy the Dirichlet boundary 
condition in the sense of viscosity solutions, as we can prove Ui satisfy the equations by 
an argument similar to the proof of Theorem 7.1. Since it is clear to see that u,i < gi on 
dfl in the classical sense from the definition of Uj, we only need to prove that (1*1,1*2) is 
a supersolution of (D e ) and particularly that u\ satisfies the boundary condition in the 
viscosity solution sense. Take xq G dfl and suppose that 

(«i-0i)(so) <0. (7.5) 

Let G C 1 (f2) satisfy {u\ — 4>)(x ) = max^W! — <p) = 0. By Proposition 7.3 we have 

ui(x ) + h 2 

. rhl\T h . 

>E i(y e~ s L^(r] h (s),-rj h (s))ds + l {h<Th} e~ h u u{h) (t] h (h)) + l {/l > T , i} e" rh 5f ly(T , i )(?7 /l (r /l ))J 
^(J^' 1 e- s L u{s) C^, -f, h (s)) ds + e-^u^hA^C^C/i A r h ))) 

for some 77^ G AC ([0, /i], fi) such that 77(0) = x and G [0, /i]. In view of (7.5), we have 
Th > for small h > 0. Therefore by a similar calculation as in the proof of Theorem 7.1 
we get 

u 1 {x )+H 1 (^,D<f>{xo)) + ^(u 1 -u 2 ){x ) >0. □ 
References 

1. S. N. Armstrong, P. E. Souganidis, Stochastic homogenization of Hamilton- J 'acobi and degenerate Bell- 
man equations in unbounded environments, J. Math. Pures Appl., to appear. 

2. S. N. Armstrong, P. E. Souganidis, Concentration phenomena for neutronic multigroup diffusion in 
random environments, submitted. 

3. G. Barles, B. Perthame, Exit time problems in optimal control and vanishing viscosity method, SIAM 
J. Control Optim. 26 (1988), no. 5, 1133-1148. 

4. F. Cagnetti, D. Gomes, H. V. Tran, Adjoint methods for obstacle problems and weakly coupled systems 
of PDE, submitted. 

5. F. Camilli, A. Cesaroni, C. Marchi, Homogenization and vanishing viscosity in fully nonlinear elliptic 
equations: rate of convergence estimates, Adv. Nonlinear Stud., 11 (2011), no. 2, 405-428. 

6. F. Camilli, O. Ley, P. Loreti, Homogenization of monotone systems of Hamilton- J acobi equations, 
ESAIM Control Optim. Calc. Var. 16 (2010), no. 1, 58-76. 

7. F. Camilli, O. Ley, P. Loreti, V. Nguyen, Large time behavior of weakly coupled systems of first- order 
Hamilton- J acobi equations, NoDEA Nonlinear Differential Equations Appl. to appear. 

8. F. Camilli, C. Marchi, Continuous dependence estimates and homogenization of quasi-monotone systems 
of fully nonlinear second order parabolic equations, preprint. 



HOMOGENIZATION OF WEAKLY COUPLED SYSTEMS OF HJ EQUATIONS 



35 



9. I. Capuzzo-Dolcctta, H. Ishii, On the rate of convergence in homogenization of Hamilton- J 'acobi equa- 
tions, Indiana Univ. Math. J. 50 (2001), no. 3, 1113-1129. 

10. G. Contreras, R. Iturriaga, G. P. Patcrnain, M. Paternain, Lagrangian graphs, minimizing measures 
and Mane's critical values, Geom. Funct. Anal. 8 (1998), no. 5, 788-809. 

11. M. C. Concordcl, Periodic homogenization of Hamilton- J acobi equations: additive eigenvalues and 
variational formula, Indiana Univ. Math. J. 45 (1996), no. 4, 1095-1117. 

12. M. C. Concordel, Periodic homogenisation of Hamilton- J acobi equations. II. Eikonal equations, Proc. 
Roy. Soc. Edinburgh Sect. A 127 (1997), no. 4, 665-689. 

13. M. H. A. Davis, Piecewise- deterministic Markov processes: a general class of nondiffusion stochastic 
models, J. Roy. Statist. Soc. Scr. B 46 (1984), no. 3, 353-388. 

14. A. Eizcnbcrg, M. Freidlin, On the Dirichlet problem for a class of second order PDE systems with 
small parameter, Stochastics Stochastics Rep., 33 (3-4):lll-148, 1990. 

15. H. Engler, S. M. Lenhart, Viscosity solutions for weakly coupled systems of Hamilton- J acobi equations, 
Proc. London Math. Soc. (3) 63 (1991), no. 1, 212-240. 

16. L. C. Evans, The perturbed test function method for viscosity solutions of nonlinear PDE, Proc. Roy. 
Soc. Edinburgh Sect. A 111 (1989), no. 3-4, 359-375. 

17. L. C. Evans, Periodic homogenisation of certain fully nonlinear partial differential equations, Proc. 
Roy. Soc. Edinburgh Sect. A 120 (1992), no. 3-4, 245-265. 

18. L. C. Evans, D. Gomes, Effective Hamiltonians and averaging for Hamiltonian dynamics. I, Arch. 
Ration. Mech. Anal. 157 (2001), no. 1, 1-33. 

19. B. Fchrman, Stochastic homogenization of monotone systems of viscous Hamilton- J acobi equations 
with convex nonlinearities, preprint. 

20. D. Gomes, A stochastic analogue of Aubry- Mather theory, Nonlinearity 15, (2002), no. 3, 581-603. 

21. K. Horie, H. Ishii, Homogenization of Hamilton-! acobi equations on domains with small scale periodic 
structure, Indiana Univ. Math. J. 47 (1998), no. 3, 1011-1058. 

22. H. Ishii, A boundary value problem of the Dirichlet type for Hamilton- J acobi equations, Ann. Scuola 
Norm. Sup. Pisa CI. Sci. (4) 16 (1989), no. 1, 105-135. 

23. H. Ishii, S. Koike, Viscosity solutions for monotone systems of second-order elliptic PDEs, Comm. 
Partial Differential Equations 16 (1991), no. 6-7, 1095-1128. 

24. H. Ishii, K. Shimano, Asymptotic analysis for a class of infinite systems of first- order PDE: nonlinear 
parabolic PDE in the singular limit, Comm. Partial Differential Equations 28 (2003), no. 1-2, 409-438. 

25. E. Kosygina, F. Rezakhanlou, S. R. S. Varadhan, Stochastic homogenization of Hamilton- J acobi- 
Bellman equations, Comm. Pure Appl. Math. 59 (10) (2006) 1489-1521. 

26. S. M. Lenhart, N. Yamada, Viscosity solutions associated with switching game for piecewise- deter- 
ministic processes, Stochastics Stochastics Rep., 38 (1): 27-47, 1992. 

27. P.-L. Lions, G. Papanicolaou, S. R. S. Varadhan, Homogenization of Hamilton- J acobi equations, un- 
published work (1987). 

28. P.-L. Lions, P. E. Souganidis, Correctors for the homogenization of Hamilton- J acobi equations in the 
stationary ergodic setting, Comm. Pure Appl. Math. 56 (10) (2003) 1501-1524. 

29. P.-L. Lions, P. E. Souganidis, Stochastic homogenization of Hamilton- J acobi and "viscous" Hamilton- 
Jacobi equations with convex nonlinearities-revisited, Commun. Math. Sci. 8 (2) (2010) 627-637. 

30. H. Mitake, H. V. Tran, Remarks on the large time behavior of viscosity solutions of quasi-monotone 
weakly coupled systems of Hamilton- J acobi equations, Asymptot. Anal., to appear. 

31. F. Rezakhanlou, J. E. Tarver, Homogenization for stochastic Hamilton- J acobi equations, Arch. Ration. 
Mech. Anal. 151 (4) (2000) 277-309. 



36 



H. MITAKE AND H. V. TRAN 



32. R. W. Schwab, Stochastic homogenization of Hamilton- J acobi equations in stationary ergodic spatio- 
temporal media, Indiana Univ. Math. J. 58 (2) (2009) 537-581. 

33. K. Shimano, Homogenization and penalization of functional first-order PDE, NoDEA Nonlinear Dif- 
ferential Equations Appl. 13 (2006), no. 1, 1-21. 

34. P. E. Souganidis, Stochastic homogenization of Hamilton- J acobi equations and some applications, 
Asymptot. Anal. 20 (1) (1999) 1-11. 

35. H. V. Tran, Adjoint methods for static Hamilton- J acobi equations, Calc. Var. Partial Differential Equa- 
tions 41 (2011), no. 3-4, 301-319. 

(H. Mitake) Department of Applied Mathematics, Faculty of Science, Fukuoka Univer- 
sity, Fukuoka 814-0180, Japan 

E-mail address: mitake@math.sci.fukuoka-u.ac.jp 

(H. V. Tran) Department of Mathematics, University of California, Berkeley, CA 94720, 
USA 

E-mail address: tvliung@matli.berkeley.edu 



